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Abstract 

Motivated by recent work on integrable flows of curves and 1+1 dimensional sigma 
models, several 0(A^)-invariant classes of hyperbolic equations Utx = f{U,Ut,Ux) 
for an iV-component vector U (f, x) are considered. In each class we find all scaling- 
homogeneous equations admitting a higher symmetry of least possible scaling weight. 
' Sigma model interpretations of these equations are presented. 

d ' 
•1—1 ■ 

1 Introduction 

. Integrability theory of nonlinear partial differential equations (PDEs) has many different 

J> I aspects which have developed over the years — Backlund transformations, Lax pairs, in- 

*^ ' verse scattering transform and soliton solutions (Calogero's "S-integrabihty" ) , linearizing 

• transformations ( "C-integrability" ) , reduction to Painleve transcendents, hierarchies of 

CN ! higher symmetries and conservation laws, master symmetries, recursion (Nijenhuis) op- 

I erators, bi-Hamiltonian structures, and others. To date the most computationally direct 

' and effective test of integrability has proved to be the condition that a PDE system should 

possess sufficiently many higher symmetries ^1 (see also [SHI E[ 1^ ) • In particu- 
lar, for all currently known examples of nonlinear scalar PDEs, the existence of one higher 
Ch ■ symmetry implies the existence of infinitely many, i.e. a symmetry hierarchy. 

^ . A rigorous proof of the observation "one symmetry implies infinitely many" was estab- 

^ I lished in [231 15] for a wide class of semilinear scalar evolutionary PDEs 

H : 

^ Ut=Unx + f{u,Ux,...,Un-lx), n > I, (1.1) 



under a homogeneity restriction with respect to a scaling symmetry on t, x, u (so-called 
A- homogeneous form). The proof uses the symbolic method of Gel'fand-Dikii, combined 
with computer algebra computations. As a main result, it was shown that any such 
symmetry-integrable scalar polynomial PDE with positive scaling weight (A > 0) for u 
belongs to one of the well known hierarchies ^ ^Jj: n = 2 Burger's (A = 1); n = 
3 Korteweg de Vries (KdV) (A = 2), modified Korteweg de Vries (mKdV) (A = 1), 
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Ibragimov-Shabat (IS) (A = 1/2); n = 5 Kaup-Kupershmidt (KK) (A = 2), Sawada- 
Kotera (SK) (A = 2), Kupershmidt (K) (A = 1), potential Kaup-Kupershmidt (A = 1), 
and potential Sawada-Kotera (A = 1). The main evolutionary PDEs of semilinear form 
not covered by this classification are the hierarchy of the nonlinear Schrodinger (NLS) 
equation and its various derivative-type generalizations in which u is complex- valued 
so its real and imaginary parts satisfy a pair of coupled scalar evolutionary PDEs. There 
has been much work and interest in finding multi-component generalizations of all these 
hierarchies, particularly where u is replaced by a vector U or scalar-vector pair {u, U) or 
matrix lA. Recent symmetry-integrability classification results in this direction have been 
obtained by extensive computer algebra computations [201 [2E1 OOj • These classifications 
test for the existence of a higher symmetry, ^ in analogy with the scalar case. 

Comparatively less work has been devoted to investigating symmetry-integrability of 
hyperbolic PDE systems jSl], apart from the sine-Gordon (SG) equation and its scalar 
variants related to the 5^ sigma model, which both have the wave equation form utx = 
f{u,ut,Ux) and possess zero scaling-weight for u. 

In this paper we present some symmetry-integrability classifications of hyperbolic vector 
PDEs. The classifications are motivated by the following three types of examples of scalar 
wave equations: 

(1) the SG equation utx = sinu, which is distinguished by scaling weight A = for 
(2a) the variant SG equation utx = U\J\ — u'f coming from the transformation u 

sin~^ (ut), which has scaling weight A = 1 for u and shares the same higher symmetries as 

the mKdV hierarchy; 

(26) the complex-valued version utx = U\/^ — {utl"^, which shares the higher symmetries 
of the NLS hierarchy; 

(3) the mKdV equation in potential form utx = u^x + u'^U2x and the NLS equation in 
potential form iutx = u^x + Wxl'^Ux- 

Our results give interesting vector counterparts of the first two types of wave equations. 
We show that these types of hyperbolic vector equations all have a natural geometric 
origin in terms of integrable nonlinear sigma models and integrable inverse flows of curves 
in Riemannian manifolds. For the third type of wave equation, we find a null result that 
no vector counterparts exist. 

2 Classification Results 

For an A^-component vector U{t,x), with N arbitrary, consider an 0(A)-invariant vector 
PDE of generalized evolutionary form 

Utmx = fiU,Ux,...,Unx,Ut,...,Utm-ix), n,m>0 (2.1) 

namely / is a sum of terms proportional to its vector arguments with scalar coefficients 
depending on dot products (•, •) of these vectors. Such a vector PDE is strictly evolutionary 

^It should be noted that there are special classes of Bakirov-type multicomponent evolutionary PDE 
systems that, surprisingly, admit only finitely many higher symmetries However, all these systems 
are of a triangular form in which one PDE is a decoupled linear evolution equation. In contrast, the 
classifications in 20 28 30 in addition to our results deal with fully coupled, nonlinear multicomponent 
PDE systems, describing a single vector or matrix nonlinear evolution equation for U or W, or a non- 
triangular scalar- vector pair of nonlinear evolution equations for (u, U). 
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if m = 0, n > 1, or is hyperbolic if m = 1, n < 1; we call it a wave equation if m = 1, n > 0. 
In all cases it is semilinear whenever n < m or if / is linear in Unx when n > m. We say 
it is of minimal differential order (m, n) provided / is not a total derivative with respect 
to X. The order of the vector PDE H2.1|) in the ordinary sense is max(n, m + 1). 

If we view a vector PDE (|2.1|) as formally defining the generator of a flow D^D^ in the 
jet space = (t, x, U, Ut, Ux, ■ ■ ■ , Uktix, ■ ■ then a higher symmetry on the solution jet 
space R^°°^ C J^°°^ (i.e. modulo equation (|2.1j) and differential consequences) will be a 
commuting flow whose generator is of the form 

Ur = g{U, Ux,..., Urx, Ut,..., Utm-lx) (2.2) 

for some order r > max(n,m + 1), where 5 is a sum of terms proportional to its vector 
arguments with scalar coefficients depending on dot products (•, •) of these vectors. Thus 
the higher symmetry condition is 

[D^, DtD^] = on . (2.3) 

A flow (|2.1j) is said to be A-homogeneous of weight (A, /i) if it admits a scaling symmetry 
group (with parameter e G M) 

{t, X, U) {e^H, e'x, e'^'U). (2.4) 

Note, with this convention for scaling weights, the order of a semilinear flow 1)2. 1|) is equal 
to the weight fi in the evolutionary case and /x + 1 in the hyperbolic case. Without loss 
of generality, higher symmetries (|2.2() for such flows can be assumed to be homogeneous 
under the same scaling symmetry group (|2.4|) (since any symmetry necessarily decomposes 
into a sum of terms with homogeneous forms which themselves will define symmetries). 
Correspondingly, we will refer to higher symmetries (|2.2|) by their weight as defined through 
the scaling 

r ^ e^'T (2.5) 

induced by the scaling symmetry group (|2.4j) . Note this weight 1/ is equal to the order r 
precisely when a higher symmetry is of semilinear form, i.e. linear in Urx. 

Throughout, we will assume a partial polynomial form for A-homogeneous vector PDEs 
H2.1|) and higher symmetries ()2.2|) . More specifically, as explained in Sec. El / and g will be 
linear in their vector arguments and polynomial in all nonnegative-weight dot products of 
these vectors, with undetermined coefficients that either are functions of any zero-weight 
dot products or otherwise are constants. We refer to such a form as the nonnegatively- 
weighted polynomial class of vector PDEs and higher symmetries. Hereafter we restrict at- 
tention to scaling weights A = 0, ^, 1, 2 in analogy with the classes of symmetry-integrable 
scalar evolutionary PDEs summarized in Sec. ^ 

In the case of evolutionary vector PDEs, the following classification was obtained in 
pSj . Up to scalings of t, x, U, the only vector counterparts of A-homogeneous evolutionary 
scalar polynomial PDEs (jl.lj) with scaling weig hts A = i, 1,2 possessing a homogeneous 
higher symmetry of least possible weight are given by, in the case of real-valued U: 

• vector mKdV equations A = 1, /i = 3, 

Ut = U3x + {U,U)Ux, (2.6) 
Ut = Usx + {U, U)Ux + {U, Ux)U; (2.7) 



4 



S.C. Anco and T. Wolf 



• vector IS equation A = ^ 




(2.8) 




(2.9) 
(2.10) 



iUt = U2cc + ai{U,U)U, + a2{U,U)a, + a3{U,U,)U + ai{U,U,)U + a5{U,U.,)U 



whose coefficients ai, . . . ,as depend on two constant parameters, falling into 6 classes, as 
listed in . We remark that this classification does not consider the complex- valued case 
for ;U = 3 or ;U = 5, i.e. vector analogs of complexly-coupled mKdV, IS, K, KK, SK, and 
potential KK, SK equations, other than ones contained in the NLS hierarchy. Hereafter, 
in the case of complexed-valued U, we allow t ^ it, r — > ir, and impose invariance under 
the phase symmetry U — > e^'^U (e G M), as displayed by NLS type systems. 

2.1 Vector-potential wave equations 

We now proceed to discuss symmetry-integrability classifications of vector Utx equations. 
Our results in theorems^lSlare established by computer algebra computations, as outlined 
in detail in Sec. 01 In each case we explicitly determine all equations admitting a higher 
symmetry of least weight. 

The first class we will consider consists of potential equations related to the integrable 
evolutionary vector equations H2.6|) - H2.11() via U = Vx, with scaling weight Ay = A — 1. 
Note that the only ones whose potential form has non-negative scaling weight Ay are the 
vector mKdV equations and NLS equations with Ay =0. A natural question is whether 
they still possess higher symmetries in potential form. We will settle a slightly more general 
classification problem by considering Ay-homogeneous vector potential wave equations 



where / is a sum of terms proportional to the vectors Vi^, < / < n, with coefficients 
depending on the dot products {Vjx,Vkx)- 

Theorem 1. In the nonnegatively-weighted polynomial class for the potential-mKdV weight 
(Ay,/i) = (0,3), every Xy -homogeneous vector wave equation i2.1S^} that possesses a ho- 
mogeneous higher symmetry of weight i/ = 5 containing no t derivatives is linear. The 
same result is true for the potential-NLS weight (Ay , /i) = (0, 2) with = 3. 

In contrast to this null result, there are many nonlinear scalar utx equations possessing 
a higher symmetry: for example, the potential-mKdV equation utx = u^x + u^Ux, and the 
3rd-order Burger's equation in potential form utx = U4x + UxU^x + ■"2^^ + \ux^U2x- 



+a6{U, U){U, U)U + miU, U){U, U)U + as{U, JjfU 



(2.11) 



Vtx = f{V,Vx,...,Vnx), n>0 



(2.12) 
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2.2 Vector sine-Gordon equations 

The next class of vector wave equations we investigate will be hyperbolic vector vari- 
ants of the SG equation utx = sinu, distinguished by the scaling symmetry {t,x,u) 
{e~^t,e^x,u). Although the transcendental form of nonlinearity in the SG equation pre- 
cludes any obvious vector analog, we note there are other integrable hyperbolic scalar 
equations with the same scaling symmetry yet compatible with a vector nature for u, 
such as the equation utx = —u{l — v?')~^utUx which is related to the scalar wave equation 
Utx = by the transformation u — > sinu. We now classify analogous hyperbolic vector 
equations possessing higher symmetries with the same weight as the ones in both the SG 
hierarchy and the hierarchy connected with the scalar wave equation. 

Theorem 2. In the nonnegatively-weighted polynomial class for the SG weight (A, fi) = 
(0,-1), every X-homogeneous nonlinear hyperbolic vector evolution equation \2.1\) that 
possesses a higher symmetry of weight u = 3 is given by (up to scalings oft,x, U) 



is an admitted higher symmetry of order 3 for weight u = 5. 

We will give a geometrical interpretation for this wave equation in terms of a A^- 
dimensional sigma model based on a paraboloid Riemannian space in Sec. ^ where we 
also discuss the question of its integrability. 

2.3 Vector hyperbolic flow equations 

The final class of vector wave equations we will consider is motivated by the SG variant 
equation utx = uy^l — whose associated flow commutes with the mKdV hierarchy of 
higher symmetry flows. With respect to this hierarchy, this wave equation represents a 
— 1 flow as it is mapped into the trivial flow Ur = hy the mKdV recursion operator. 
We recall that the recursion operator generates the mKdV hierarchy starting from the 
j;-translation symmetry Ur = Ux, referred to as the flow for the hierarchy; the +1 flow 
is the mKdV equation itself Ur = u^x + ^u'^Ux, while its higher symmetries define the +2 
flow and so on. We now systematically look for hyperbolic vector equations that describe 
— 1 flows for the hierarchies of each integrable vector evolutionary equation ()2.6() - (|2.11() . 
In general, the —1 flow corresponding to a A-homogeneous +1 flow in the vector case will 
have weight (A, —1), as we explain next. 

Consider first the case of the mKdV and IS hierarchies where the flow is Ur = Ux 
with scaling weight vq equal to 1. Since the scaling weight of the +1 flow is z^i = 3, the 
weight difference between successive flows in the hierarchy is w = z^i — i^o = 2. Hence if 
the —1 flow is regarded as a nonlocal evolutionary equation Ut = D^^ f{U,Ut,Ux) then 
it will have scaling weight ^ = vq — w = 2vq — vi = —1. The remaining case is the 
NLS and derivative-NLS hierarchies where the flow is Ur = iU with scaling weight 




a = -1,0,1 



(2.13) 
(2.14) 



This result extends at least to weights = 4, 5, 6, 7, and in particular 



Ur = m Usx) + S{Ux, U2x)){U, Ux)Ux + (([/, U2x) + {Ux, Ux)fUx 
+a{{Ux,U3x) + {U2x, U2x))Ux 



(2.15) 
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uq equal to 0. The +1 flow is then Ur = Ux, and the (derivative-) NLS equation itself 
corresponds to the +2 flow, with scaling weight z/2 = 2. The weight difference between 
flows in this case is = i^i — fo = 1, and so the —1 flow will then have scaling weight 
= — w = 2vq — vi = —1. Therefore in all cases the weight // = —1 of t for the —1 
flow is simply the negative of the weight of x. A similar argument clearly applies to the 
well-known scalar hierarchies. 

We now state classifications that determine all —1 flows in the vector case. 

Theorem 3. (rriKdV and NLS —1 flows) For X-homogeneous nonlinear hyperbolic vector 
evolution equations \2. 1]) with weight (A,^) = (1,-1) in the nonnegatively-weighted poly- 
nomial class, the ones that possess a higher symmetry of mKdV weight u = 3 are given by 
(up to scalings oft,x,U) 



Utx = ±^l-a{Ut,Ut)U, a = -1, 1 (2.16) 
Ur = U:ix + \oi{U,U)Ux, (2.17) 



and 



Utx = {U,Ut){UuUt)-\l + A)Ut-\a\{l-A)U, a = -1,0,1 (2.18) 

Ur = [/3x + 3a(f/, U)Ux + 3a{U, Ux)U, (2.19) 

where A = ±^/l^^^^^(JJuU^ , (2.20) 

while the ones that possess a higher symmetry of NLS weight u = 2 are given by ( up to 
scalings of t,x,U) 

Utx = {U,Ut){Ut,Uty\a + B)Ut-{a-B)U, a = -1,1 (2.21) 

iUr = U2x + ^{U,U)U (2.22) 



where B = ±\j - ^{Ut, Ut), (2.23) 



and 



^Jc2 + {Ut,Ut) + a Cf 1 {{U, Ut)Ut - {U, Ut)Ut) + C2 - {UuUt) - a U, (2.24) 
iUr = U2x + l{U,U)U -^{U,U)U, a = -1,0,1 (2.25) 



where Ci = J {Ut,Ut){Ut,Ut) - {UuUt)\ (2.26) 



C2 = ±a/«' + MUu Ut) + {Ut, Ut){Ut, U 



th 



All four of these —1 flow equations are expected to be integrable, as they possess a 
-|-1 flow of the vector mKdV hierarchies or a +2 flow of the vector NLS hierarchies (up 
to scaling of t, x, U,t). We will interpret the flrst two of them geometrically in terms of 
homogeneous N + 1-dimensional gauged sigma models in Sec. IH 

Finally, we mention that a classiflcation of IS and derivative-NLS —1 flows analogous 
to theorem|21in the nonnegatively-weighted polynomial class with weight (A, /i) = (^, —1) 
yields a null result. However, such flows do exist in a wider class, which we will present 
elsewhere. 



Some symmetry classifications of hyperbolic vector evolution equations 



7 



3 Computational Aspects 

Here we describe the algorithmic and implementation issues involved in the computations 
used to prove the classification theorems in Sec. IS These computations divide into three 
parts: 

(1) making an ansatz for the PDEs and symmetries with specified homogeneity weights; 

(2) computing the symmetry conditions; 

(3) solving the system of symmetry conditions. 

The first two tasks are handled by the program SVSym (Scalar- Vector-Symmetry) ^ and 
the third computation is performed by the package Crack described in j32j . 

The class of PDE systems for investigation is allowed to involve any number of scalar 
functions u and vector functions U of t,x, which we denote collectively as (/>*. An essential 
limitation currently in the program is that the left-hand side (l.h.s.) of each PDE must 
consist of a single derivative (pptq^ sharing the same order p,q > for all PDEs in the 
system. 

It is assumed that all derivatives occurring on the right-hand side (r.h.s.) of the PDEs 
and the symmetries are not equal to any l.h.s. derivatives or their differential consequences. 
In addition, r.h.s. derivatives must have a lower priority with respect to some total or- 
dering >T of derivatives than all l.h.s. derivatives, to avoid infinite substitution loops in 
computing the symmetry condition. For example, the equation utx = uu + Uxx would 
lead to an infinite loop of substitutions for differential consequences of utx, as seen by 
utxx — *■ uttx + Uxxx, uttx uttt + utxx- This problem is avoided if we adopt (without 
loss of generality) the lexicographical ordering t >t x. Then because we have uu >t utx, 
this equation would be handled by bringing the highest priority derivatives to the l.h.s, 
utt = Uxx — Utx, which no longer leads to any infinite substitution loops. 

3.1 Specifying an Ansatz 

Many integrable systems are A-homogeneous with respect to multiple homogeneity weights. 
For example, denoting the weight for x as ^, the vector NLS system |2H1 

Ut = Uxx + {U, V)U, Vt = -Vxx - {U, V)V 

and its higher symmetry hierarchy is A-homogeneous with weights /i,^, Ac/,Ay equal to 
2, 1, 1, 1, as well as 0, 0, 1, —1. The latter weights have been used as an extra filter in our 
investigation of PDEs in the complex- valued case. To accommodate multiple weightings, 
the program SVSym allows the specification of a list of weight sets ^, Aj}. 

The program places two a priori restrictions on possible weights to ensure that the 
number of possible homogeneous terms allowed to appear in the PDEs and symmetries is 
necessarily finite. A first restriction is that x cannot have zero weight in all weight sets 
as this would permit an unlimited differential order of x- derivatives (note the differential 
order of t- derivatives is always finite due to the ordering t >t x combined with the fixed 
number p of t- derivatives on the l.h.s.). A second restriction is that there must be a 
weight set {fi* , , u* , X*} for which the total weight of the l.h.s. of the PDEs (pptgx and 
symmetries 4>l. is non-negative, namely wf^^ := A* + pjJ* + qC > (with 7^ 0), and 

^ SVSym is a completely new program, more general than the earlier program used for the classification 
of evolutionary PDEs in 1281 13U| . which was limited to positive weights and a strictly polynomial ansatz. 



8 



S.C. Anco and T. Wolf 



y^sym ^* _^ ^* y then the r.h.s. of the PDEs and symmetries can be generated by 
the following finite ansatz: 

<,. = r, U;,^, = Y,F^''UI and = g\ ^ = Y,G^''u'j. (3.1) 

k,J k,J 

Here J is a jet index, denoting partial derivatives including those of zeroth order (only t- 
and x-derivatives of u^, are used whose total ordering with respect to t >t x is lower 
than that of the l.h.s. of the PDEs). In (jXTJ), [/) has *-weight > 0; f\g\F'^-^,G'''-^ are 
the most general polynomials (of appropriate weight) that can be built from the following 
scalar polynomial variables, all of which are constructed to have positive *-weight: Uj, 
products Uj^Uj^ where one of the two factors has a negative *-weight, and dot products 

{Uj^,Uj^) where at most one of the two factors may have a negative *-weight. The 
undetermined coefficients in all polynomials are arbitrary functions of similarly constructed 
scalar variables Uj, Uj^u^^^, {Uj^,Uj^) whose weight is zero in all weight sets. All weight 
sets other than the *-weight set are used as an extra filter on the generated terms. (The 
restriction of using polynomial variables built from at most two scalar or vector factors 
could obviously be generalized to allow more factors such that the total weight of the 
product still satisfies the previous conditions.) 

It is important to emphasize that the ansatz (|3.H) for the PDEs and symmetries need 
not be strictly polynomial if /i* < or ^* < or any X* < 0. 

A few flags give convenient means to cut down the generality of the ansatz, for example 
by discarding t-derivatives or imposing constant coefficients in the PDEs or symmetries 
or both. Extra conditions or inequalities on coefficients can be added easily. 

The following examples show the ansatz generated for the classifications in theorems I^J 
and|31 The vector SG ansatz with fi = —1, u = 3, Xjj = is: 

Utx = aiU, 

Ur = hUsx + b2U2AU., U) + b3U^{U2^, U) + hdJ^Ux, U^) + hU^{U, + beU{Us^, U) 
+hU{U2^, U^) + bsU{U„ U^){U, U,) + bgU{U2^, U){U, U,) + bioU{U,, uf 

where all coefficients Oj, bj are undetermined functions of {U, U), {Ux, Ut). The ansatz for 
the vector mKdV inverse flow with ^ = —1, v = 3, Xjj = 1 is: 

Utx = aiUt{U,Ut)+a2U, 

Ur = biU3x + b2U2AU, Ut) + bsUxiUx, Ut) + b^UxiU, Ut)^ + b^UxiU, U) + beUtiUx, Ut)^ 
+bjUt{Ux, Ut){U, Ut)^ + bsUtiUx, Ut){U, U) + bgUtiU, Ut)\U, U) + bioUtiU, Ut)^ 
+biiUt{U, Ut){U, Ux) + buUtiU, U2x) + bi3Ut{Ux, Ux) + b^UtiU, U)^ 
+bi5U{Ux, Ut){U, Ut) + bieU{U, Ut){U, U) + b^UiU, Ut)^ + bisU{U, Ux) 

where all coefficients Oj, bj are undetermined functions of {Ut, Ut). The ansatz for the vector 
NLS inverse flow is similar, using a pair of vectors U, V, with fi = —2, i/ = 6, A[/ = Ay = 1, 
and imposing the additional the filter fj, = ly = 0, Xu = 1, Xy = —1. 

3.2 Computing Symmetry Conditions 

In the special case of evolutionary PDEs, all substitutions of (pl and (/>^ in the symmetry 
conditions (/>n , := D^cpl — D^(j)\. = are done only once because the r.h.s. of the PDEs and 
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symmetries do not contain t- or r- derivatives. As a consequence, the symmetry conditions 
are separately linear in the undetermined coefficients of the PDEs and of the symmetries. 

The situation is different when the l.h.s. of the PDEs contain an x-derivative of 0*, 
such as in the hyperbolic case cpl^. Then any t- or x-derivatives that appear in the r.h.s. of 
the PDEs and symmetries lead to (repeated) substitutions of (f)l^ through the PDEs. The 
symmetry conditions become polynomially non-linear in the coefficients of the PDEs but 
stay linear in the coefficients of the symmetries. Due to this non-linearity, which is fur- 
ther amplified whenever the coefficients are functions of f-derivatives of 0* , the symmetry 
conditions can become extremely large. The size of intermediate expressions can become 
exceedingly high during the repeated substitution and simplification process. To get over 
this memory hurdle, the linearity in the coefficients of the symmetry is exploited. If the 
r.h.s. of the symmetry consists of s terms with undetermined coefficients bi, . . . ,bs then 
the full symmetry condition is computed via 



where g'^^^ denotes the kth term in the r.h.s. of (pi = , i.e. all coefficients except bk are 
replaced by zero. This trade-off of computing memory for computing time is also justified 
by the fact that the time to simplify large expressions grows non-linearly with their size, 
at least in the computer algebra system Reduce that is used. 

3.3 Solving the Symmetry Conditions 

Compared to the case of evolutionary PDEs, the symmetry conditions for the general case 
of PDEs with l.h.s. (pptqx purely algebraic but may involve ordinary and partial 

differential equations and may be inhomogeneous and non-linear in the undetermined coef- 
ficients of the PDEs. The system of symmetry conditions is fed into the computer program 
Crack, which was used to solve the algebraic symmetry conditions in the classification 
of evolutionary vector equations |2H| and coupled evolutionary scalar- vector systems [30] . 
Originally, Crack has been developed as a solver for overdetermined PDE systems, and 
this functionality is essentially required for solving the polynomially nonlinear differen- 
tial system of symmetry conditions that occur in the classifications of vector hyperbolic 
equations in Sec. [21 Solutions of non-polynomial form, like square roots in the equations 
found in theorem 121 may appear and often are in reach of the Reduce ODE solver that 
is employed by Crack when appropriate. 

The strength of Crack compared with other computer algebra packages for solving 
overdetermined systems lies in the flexibility of its approach, the possibility to run it 
in a fully automatic or interactive mode, and its support in handling especially large 
systems (safety, process control, heuristic algorithms to cut down the size of equations, 
ability to take advantage of the linearity of coefficient functions of the symmetries). In 
addition. Crack provides a verbose mode reporting all computational steps, which allows 
in principle the solutions to be checked by (human) inspection (in practise this is hard 
due to the high number of steps). 



s 




(3.2) 
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4 Sigma Model Interpretations 

The wave equation utx = and SG equation utx = smu are well-known to be intimately 
related to two of the simplest sigma models 21 , which are described by a geometrical 
wave equation for maps on 2-dimensional Minkowski space M^'^ into target spaces and 
S^, respectively. In general a sigma model on M^'^ is given by the nonlinear wave equation 



= ^v,dy = d,dy + T^^ci^)d,u^d,u^ (4.1) 

for scalar functions u'^{t, x) representing a map into some Riemannian target space (M, g^^^), 
A = l,..., dimM. Here ^V^ = d^ + T'^g^{u)^^u^ is the pullback to M^'^ of the (torsion- 
free) covariant derivative determined by the metric in local coordinates on the 
manifold M, given in terms of the Christoffel symbol ^^bc metric. This wave equa- 

tion has the geometrical action principle S[u] = J {d^u, d^u)gdtdx where (•, ■)g denotes the 
inner product in T{M) given by the metric g^^i''^)- Another geometrical meaning for this 
wave equation arises if we view u^{t, x) for fixed f as a curve embedded into the manifold 
M. The sigma model thereby describes a flow of this curve, whose arclength is preserved 
due to the conservation law 

d,{d^u^dyg^siu)) = 0. (4.2) 

Because of the reflection symmetry t ^ x in the model, there is also a conservation law 

9,(a,n^a,n%5(n))=0. (4.3) 

These conservation laws are connected with invariance of the model under conformal 
scaling transformations t ^ t' := a{t), x — > x' := P{x). This freedom can be fixed in a 
natural manner by putting a'(t)^ = {d^u,d^u)g and /3'(x)^ = {d^u,d^u)g. Finally, note 
that in general the sigma model equation 1)4. 1() is A-homogeneous with scaling weight A = 
on while t and x have opposite weights of = —1 and +1, in the conventions of Sec. IS 

0(iV)-invariant vector wave equations arise from sigma models in various ways and for 
certain types of target spaces. One direct way is by considering target spaces embedded 
as an 0{N) symmetric hypersurface in Euclidean space M^"^^ or Minkowski space R^'^. 
Vector coordinates for the hypersurface given by its embedding lead to natural vector 
generalizations of the variant scalar wave equation utx = —u{l — v?')''^utUx which is the 
sigma model for 5^ C M^. A second way involves introducing a flat Cartan connection and 
frame bundle on homogeneous target spaces with the tangent space structure g/f) ~ M^+^ 
based on pairs of Lie algebras g D [). For a suitable geometric choice of frame, the 
connection components yield natural vector generalizations of the variant SG equation 
Utx = u^Jl — uf which comes from the sigma model for 5^ ~ S0{3)/ SO{2). 

In Sec. 14.11 we show that the vector wave equation (|2.13|) in theorem |21 arises from 
the nonlinear sigma model for an 0{N) symmetric paraboloid hypersurface in M^"*"^. In 
Sec. 14.21 we derive the first two hyperbolic vector equations in theorem |31 from the flat 
Cartan connection of gauged sigma models viewed as a flow of curves on homogeneous 
spaces (Klein geometries) associated with the Lie algebras so(A^+l) C so(A^-|-2),sii(A^+l). 
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4.1 A^-sphere and A^-dimensional paraboloid sigma models 

Consider firstly the A^-sphere target space as an embedded hypersurface in M^"^"^. We 
describe by the A^ + 1-component unit vector u\ i = 1, . . . , A^ + 1, constrained to 
satisfy 1 = u^u^5^j. If we resolve the constraint by expressing u^~^^ = ity^l — w^u^dj^^ 
with = {u^,u^~^^), A = 1, . . . ,N, then represents a set of coordinates on . In 
particular the upper or lower hemisphere of is mapped by to the open unit ball 
in M^. In these coordinates, the Riemannian metric on as induced by the Euclidean 
metric 5^^ on M^^"-*^ is given by 

9ab («) = ^AB + (1 - u'^Uc)~^Uj^Ub (4.4) 

with the inverse metric g^^ = — u^u^ , where Uj^ := dj^^u^. The Christoffel symbol 
and covariant derivative on are readily found from the condition f^g"^^ {u) = 0, 
giving 



Hence the curvature is given by 

'yc[A 



"V^/V^] := i?^^^^ = 2gc.^{n)6^.^x, X = 1 (4.6) 



where x denotes the scalar curvature. Note this indicates is an 0{N) symmetric space 
of (positive) constant curvature. 

The sigma model equation (|4.1|) for the target space 5"^ thus has the form 

= d^d^u-^ + u^id^u^d^uB + (1 - u^usy^ucu^d^u^d^u^) (4.7) 

which is an 0(A^)-invariant hyperbolic vector equation for u^{t, x). The scalar case A^ = 1 
is seen to reduce to = utx + w(l — v?)~^utUx- Integrability of this sigma model for all 
A^ > 1 is a well-established result via a Lax pair associated with 5''^ = 5'0(A^ + 1)/50(A") 
as a Riemannian symmetric space |33| I21j , and also more recently via a bi-Hamiltonian 
structure (higher-symmetry recursion operator) derived from a moving frame formulation 
of as a constant curvature space |24| |2j. Its higher symmetries were first found in 
j21j using a Backlund transformation technique which relies on conformally scaling t — > t', 
X ^ x' through the conservation laws (|4.2j) and (|4.3)1 . The hierarchy of higher symmetries 
is homogeneous and has positive scaling weight u = 3, 5, . . ., with respect to the scaling 
group {x' ,t\u^) — > {e'^x' ,e~''t' ,u^), e G M. But when expressed in the original variables 
i, X, the higher symmetries are all homogeneous of weight i/ = 0. For this reason the vector 
wave equation ()4.7() does not appear in our classifications. 

To account for the vector wave equation (|2.13|) and its higher symmetry (|2.15|) with SG 
scaling weight in theorem [21 we look for a target space with an appropriate Riemannian 
metric. Working backward from the nonlinear terms in the equation 

= d^d^u"^ + u^(l + u^ucy^^^u^^.^UB (4.8) 

we deduce the Christoffel symbol 

T^Bci'^) = u^Bcil + u^uj^r^ (4.9) 
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and hence 



= + ""^^Bci^ + ^^^d)-' (4.10) 



is the covariant derivative. Prom the condition ^g^^ (u) = 0, it is straightforward to 
derive g^^ = 6^^ + (1 + Uq)~^u'^u^ which then gives the metric 

9ab (^) = ^AB - UaUb- (4.11) 

Finally we are able to identify this metric as the induced Riemannian metric on the A^- 
dimensional target space given by the paraboloid hypersurface 

^N+i ^ In^^fij^^ in M^+i. (4.12) 

This target space is an 0{N) symmetric Riemannian manifold, whose scalar curvature is 
non-constant 

X = (l + 2u^+i)-i, R^sc'' = '^^C[A9B]''i^)x (4.13) 

where gj^^ := dj^fjg'^^. (Curiously, the metric tensor gj^^{u) on the A^-dimensional 
paraboloid is the same as the inverse metric tensor g^ ^{u) := g^'^ {u)6(jj^ on the N- 
sphere.) 

We note there is a Lorentzian version of these sigma models 1)4. 7p and ()4.8() obtained 
via replacing the Euclidean space R-'^"'"^ by Minkowski space M'^'^. Correspondingly, 
is replaced by its Lorentzian counterpart given by the A^-dimensional hyperboloid space 
u^~^^ = iby^l + w^u^S^^ which comes from the spacelike hypersurface 1 = u^u^r]- in 
R^'ij where rj-j is the Minkowski metric and = {u'^,u^~^^). Its induced Riemannian 

metric is g^si'^) ~ ^AB + (-*- + ^'^^c)~^'"a^B' with constant (negative) scalar curvature 
X = —1. In contrast, the iV-dimensional paraboloid space u^~^^ = ^u^u^6^^ is its own 
Lorentzian counterpart, but with the metric replaced by g^^ (u) = 6^^ +u^u^ as induced 
by the Minkowski metric rj-, and with non-constant scalar curvature X = (1 ~ 2u^~^^)~^ . 
Note in the Lorentzian case this paraboloid hypersurface has a null tangent direction at 
points 2u^^'^ = 1 in M.^'^ where both the metric and curvature become singular. 

Thus the nonlinear sigma model 1)4. 8() for the A^-dimensional paraboloid target space 
in the Euclidean and Lorentzian cases reproduces the vector wave equation (|2.13|) with 
the vector U identified as the local coordinates u"^. Integrability of this model does not 
seem to be known in the literature, so its higher symmetry 1)2. 15|) stated in theorem |21 is 
a new result. Unlike the model, the existence of a hierarchy of higher symmetries is 
an open question, apart from the scalar case = 1 of this model ()4.8|) which gives the 
same integrable sigma model as the scalar reduction of the model (|4.7|1 . As noted 
in the theorem, for A > 1 we checked that there are higher symmetries of all weights 
u = 4,5, 6, 7. This strongly suggests the model (|4.8j) is indeed integrable. A proof will 
require exhibiting e.g. a Lax pair, a bi-Hamiltonian structure, or a symmetry-recursion 
operator. 
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4.2 Homogeneous (gauged) sigma models and curve flows 

We now turn to the homogeneous target spaces M = G/H based on Lie groups H C G 
chosen to be semi-simple and compact, with the structure of a symmetric space. This 
structure is expressed by the Lie algebra decomposition g = f) © p such that [f), p] C p and 
[pjp] ^ f), where \},q denote the respective Lie algebra of H,G, and p is a vector space 
canonically identified (through left-invariant vector fields) with the tangent space TM at 
any point in M. As is well-known, these target spaces carry the structure of a symmetric 
Riemannian manifold whose metric tensor is given by the Cartan-Killing inner product 
(•, •) on TG ~ Q restricted to TM ~ p. This metric and its curvature tensor are both 
G-invariant and covariantly constant. A few main examples are SO{k + 1)/S0{k) ~ S^, 
SU{k + l)/U{k) ~ CP'', SU{k)/SO{k) (see dB] for a complete list and properties). The 
sigma model 1)4. If) for all such target spaces is a gauged chiral equation [5j 

= d^iu-^d^u)"^ = ^V^d^u^ (4.14) 

with u'^{t,x) given by gauge equivalence classes of G- valued functions under the gauge 
group H C G, where denotes the gauge covariant derivative. Viewed as a flow of curves 
in M = G/H, this model has a natural frame bundle formulation employing a flat Cartan 
connection (identified with the Maurer-Cartan form on G) whose components will satisfy 
an 0(A^)-invariant nonlinear vector wave equation j7j. Details of Cartan connections and 
frame bundles for homogeneous spaces are presented in j25j . 

Firstly, we consider TM ~ R^+i = p as given by g = so{N -\- 2) and f) = 5o{N + 1). 
Recall, so{k) is a real vector space of dimension ^^^^"^^ isomorphic to the Lie algebra of 
kxk skew-symmetric matrices. So TM = so{N + 2)/so{N + l) is isomorphic to p = 
as described by the following canonical decomposition 



-p h 



°T l)^so{N + 2), hGso(iV + l), pG]R^+^ (4.15) 

A 



parameterized by the N + 1-component vector p , A = 1, . . . , N + 1. Note the Cartan- 
Killing inner product on q is (up to a normalization factor) the trace of the product of 
5o{N + 2) matrices 1)4. 15() and hence reduces for h = to the ordinary dot product 6^^ on 
vectors p'^. On the target space M = SO{N + 2)/ SO{N + 1) we introduce a flat Cartan 
connection ^lij^ which is so{N + 2)-valued with the decomposition ^o;^ = ^cOj^ + ''i^^- 
As a consequence of the Lie algebra structure ()4.15|1 . ^lo^ provides an 5o{N + l)-valued 
connection to which is associated an M^^^-valued orthonormal frame e for the tangent 
space TM ~ 5o{N + 2)/so{N + 1) ~ M^+^. More precisely, the underlying (gauge) 
structure group here is SO{N + 1), namely the matrix Lie group associated with i) = 
5o{N + 1). Up to the gauge freedom of SO{N + 1) rotations, there is a canonical soldering 
relation fu;^ = between the frame and the reduced connection, where is the coframe 
satisfying (e^, e^)p = . This soldering can be shown to encode a nonzero Riemannian 
curvature of the metric on M, determined by the coframe and the connection ^w^- 

To derive an 0(A^)-invariant vector wave equation from the gauged sigma model (|4.14|1 
for this target space M = SO{N + 2)/S0{N + 1), we first write ^uot ■= ^uj^^tU^, ^^x ■= 
^ujj^d^u^, et := Cj^d^w^, Cx ■= Cj^d^w^, in terms of u"^{t,x). These obey 

= d^Sx + ^uJtex = d^et + ^w^e* (4.16) 
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as a consequence of the chiral equation (|4.14j) and the structure equation of the soldering 
frame = e^. Interpreted geometrically, the second equality is the Cartan structure 
equation for vanishing torsion of the connection ^w^. The Cartan structure equation for 
the curvature of this flat connection gives 

djiot - d^'^LO^ + i'^iu^.M] = = et A ex (4.17) 

where A denotes the antisymmetric tensor product. Note from equation (|4.16jl we recover 
the conservation laws 1)4. 2|) and ()4.3|) for this sigma model, = d^{ex,ex)p = d^{et,et)p. 
By conformally scaling t,x, we normalize {ex,ex)p = {et,et)p = 1. Then we use the 
inherent gauge freedom on using the rotation group SO{N + 1) to put 

6x^+^ = 1, 6x^ = 0, (4.18) 

without loss of generality. Thus one of the conservation laws is satisfied identically; the 
other one determines 



et''^' = ±^Jl-et'e/6^^. (4.19) 

In the index notation here, Cx = (e^*, Cx^"*"^) and et = {et^ ,et^~^^), i = l,...,N, are 
(iV + 1 X 1) column matrices comprising the frame components of, respectively, the tangent 
vector d^u^ along the curve u^{t,x) and the direction vector d^w^ of the flow on M = 
SO{N+2)/ S0{N+1). The condition ()4.18|) geometrically means restricted to the curve 
u^{t,x) is an adapted moving frame, consisting of N normal vectors and one tangential 
vector (so x has the geometrical meaning of arclength). There now remains SO{N) gauge 
freedom that preserves the adapted form (|4.18j) of ex, which is described by rotations 
in the normal space along the curve. We use this freedom to fix the connection ^u!x 
so it corresponds geometrically to specializing the adapted moving frame to be parallel 
(cf. 6 ) with respect to the curve u^{t,x) as follows. Write ^uJx = (''wx*'', ^'^x*^'*'^) and 
^ujt = i^ijJt\ ^ujt^^^) for the components of the (A^ + 1x^ + 1) connection matrices. 
The curvature equation ()4.17|) shows that the pair of matrices ^lVx^'' , ^'^t'' constitutes a flat 
so(A^)-valued connection 1-form, and consequently along the curve we put ^Wx*"* = by the 
available gauge freedom. This leads to ^uJt'' = without loss of generality. Geometrically 
speaking, ^lOx then describes the derivative of e under infinitesimal displacement along 
the curve u^{t^x) such that the normal vectors in the frame have a purely tangential 
derivative. 

Now, the torsion equation (|4.1()j) can be used to find 

fa;/^+i = 0, (4.20) 
^u;x^^+^ = (e,^+^)-ia,e/, (4.21) 

while from the curvature equation (|4.17|1 we have 

e/ = -d^'^u^x'"'^^. (4.22) 

Thus the reduced Cartan equations (|4.21|1 and (|4.22j) together with the relation (|4.19|) yield 
a system of coupled vector equations d^et = ±(1 — etet''5^Y/'^^uJx^^~^^ and d^^LOx^^~^^ = 
—et . Viewing ^Wx*^"*"^ as an A^-component vector U{t,x), and eliminating e^* from these 
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equations, we obtain precisely the 0(A^)-invariant nonlinear vector wave equation (|2.1(i|l 
in theorem |3J As shown by the theorem, this vector wave equation possesses a higher 
symmetry given by the +1 flow in the vector mKdV hierarchy (|2.H|) . 

Integrability of the nonlinear vector wave equation H2.16() was first shown via a Lax 
pair ini El- Related derivations of it using both a zero-curvature and a Lie-algebraic 
formulation of the Lax pair were obtained in j221 [71 131| . Our derivation here is more 
intrinsically geometrical, based on the approach in \17\ \24\ [2j relying on use of a parallel 
moving frame and a corresponding adapted connection 



e,. = (0,...,0,1) GM^^+\ V = (^_^T oJeso(iV+l), U eR'' (4.23) 

for a curve u^{t, x) associated to the sigma model (|4.14|) as interpreted as a flow of curves 
on the Riemannian space SO{N + 2)/ SO{N + 1) ~ S^~^^. The connection is related to 
the frame by the Frenet-Serret formula ^V^^e^ = — ^w^e^. In particular, in this frame the 
connection components U geometrically represent differential invariants of the curve. 

The second 0(A'^)-invariant vector wave equation (|2.18j) in theorem 01 has a similar 
derivation, starting from the target space M = SU{N + 1)/S0{N + 1). One notable 
difference in this sigma model is that the target space has dimension greater than A^+1 and 
hence a reduction condition must be imposed on the Cartan connection and its associated 
soldering frame, allowing them to be expressed in terms of a + 1-component vector. The 
higher symmetry structure for this vector wave equation obtained by theorem |31 identified 
as the +1 flow in the vector mKdV hierarchy (|2.7j) . is a new result. Of course, a proof of 
its integrability will require e.g. checking that it admits the symmetry-recursion operator 
of the hierarchy. The same remarks apply to the other two vector wave equations (|2.21|) 
and fn^i in theorem El 

To outline the derivation of equation 1)2. 18() . recall su{k) is a complex vector space 
isomorphic to the Lie algebra oi k x k skew-hermitian matrices. The real and imaginary 
parts of these matrices belong to the real vector space so(fc) of skew-symmetric matrices 
and the real vector space 5{k) ~ 5u{k) /5o{k) defined hy kxk symmetric trace- free matrices. 
Hence g = su{N+l) has the decomposition g = f)+ip where f) = so(iV+l) and p = s(A^+l). 
There is a natural subset of matrices in 5u{N + 1) given by the special form 

^ p-pI-p^pe5{N + 1), (4.24) 



N + 1 



which is parametrized by the arbitrary iV + 1-component vector p'^. This set of symmetric 
trace-free matrices ()4.24|) . hereafter denoted 5p{N + 1), provides a canonical reduction 
of 5u{N + l)/5o{N + 1) to M^"*"^ as a + 1-dimensional manifold. The restriction of 
the Cartan-Killing metric of 5u{N + 1) to 5p{N + 1) yields a Riemannian metric (•, ■)sp 
proportional to the trace of the product of matrices (|4.24|) . In particular, the norm of 
a matrix (|4.24|) in this metric (with a suitable normalization) is simply p ■ p. (Note 
the nonlinearity of the matrix norm in terms of p is due to the lack of a vector space 
structure for 5p{N + 1).) We now introduce an su{N + l)-valued flat Cartan connection 
0cj^ = ^LO^+i^iOy^, with ^u;^ defining an 5o(A^+l)-valued connection to which is associated 
an 5p{N + l)-valued soldering frame and a corresponding coframe = such that 
(e^,e^)sp = S^^. There is SO{N + 1) gauge freedom inherent in and ^uj^. This 
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structure is a reduction of the full soldering frame and flat Cartan connection attached to 
the target space M = SU{N + l)/SO{N + 1). 

To proceed we impose this reduction on the target space in the gauged sigma model 
1)4. 14(1 for u'^{t,x). Thus, along the tangent vector d^u"^ and the flow direction vector 
d^u"^ for the curve u'^it, x), we have the frame matrices e^, et and the connection matrices 
^iOx , ^ujt that obey the torsion and curvature equations 

= df^Bx + ^iOtex = d^et + ^ujxet, (4.25) 
dju;t - d^u^x + [''c^x.^u;*] = rc^x/c^t] = \ex.et\- (4.26) 

Using the sigma model conservation laws (|4.2j) and (|4.3j) . we put (e^;, ex)s,j, = (e* , ^t)^^ = 1- 
We further fix ex-.'^'-j-^x (via the SO{N + 1) gauge freedom) so they represent a parallel 
moving frame and an adapted connection along the curve u^{t,x). Namely, in matrix 
notation. 



^ G 5p(iV + l), "u;, = (^_° T o) G so(iV+l), C/ G (4.27) 



iV + 1 \0 
where S 5o{N), and also 

= TTTT^ - (Jt ] e ^^(iV + 1), {V,v) G R''^' (4.28) 



N + l v-'V'V, 
where 

v'^ + y.V = v. (4.29) 

Note 

Cx 1 s-re of the form H4.24|) with respective parameterizations chosen in vector index 
notation as 1 = p^"*"^, = p*, and y/v = p^"*"^, {\/v)~^V^ = p"; the norm of et is 1 = 
(p^"*"^)^ +p^p^S^j = V + v~^V^V-' d^j. Now, from the 5o{N) part of the curvature equation 

(|4.26|) together with the first half of the torsion equation (|4.25|) . we find ^uJt ~ q 

Then the remaining parts of equations H4.25() and 1)4. 26^ reduce to 

d^V = vU -y-^U -V V, dfU = -V, (4.30) 
d^v = -2U-V, d^{v-W^V) =V"^U + U^V. (4.31) 

We see equations (|4.30|) and (|4.29j) are equivalent to the 0(A^)-invariant nonlinear vector 
wave equation (|2.18|) . while the last equation ()4.31|) is found to hold as an identity. 



4.3 Other models and concluding remarks 

The sigma models connected with the two 0(A^)-invariant nonlinear vector wave equations 
(|2.16|) and (|2.18|) in theorem |31 exhaust all the types of real symmetric Riemannian spaces 
G/SO{N + 1). The other two vector wave equations ()2.21|) and (|2.24|) . which in contrast 
are complex-valued and C/(A^)-invariant, might then be expected to arise from hermitian 
symmetric spaces G/U{N + 1). However, sigma models based on the latter type of target 
space naturally take the form of coupled scalar-vector wave equations, and hence degener- 
ate cases would have to be investigated in order to account for these purely vectorial wave 
equations (|2.21|) and (|2.24j) . 
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Interestingly, the last vector wave equation (|2.24|) given in theorem |S1 stands out from 
the others in one respect; it is neither well-defined in the scalar case = 1 nor in the 
real- valued vector case (precisely when the expression (|2.2H|) vanishes). Thus it exhibits 
an intrinsic complex-valued vectorial nature. Some insight on this situation comes from 
the NLS higher symmetry p. 25)1 admitted by this vector wave equation. In |1^ the 
corresponding vector NLS hierarchy 1)2. 1U() is shown to originate from the matrix NLS 
equation idJJ = d'^U + UU'^U where the matrix U = u^^a is chosen to belong to the span 
of the gamma matrices 7yi, A = 1, . . . ,N, defined to satisfy the Clifford algebra relations 
7A7B + IB^A = 25^^I. This observation suggests the vector wave equation p.24j) may 
arise in a similar fashion from a suitable matrix type of sigma model. 

Finally, in connection with these remarks, we mention the deep work j27| I12L of 
Svinolupov and Sokolov on the construction of integrable vector and matrix equations of 
hyperbolic and evolutionary type, related to special kinds of nonassociative algebras. 
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